\begin{verbatim}
> p := 2^131+1883;
> fac := Factorization(p-1);
> TrialDivision(p-1);
[ <2, 1>, <3, 1>, <5, 1>, <37, 1> ]
[ 2452485527358115051988285458967698823 ]
> Factorization(2452485527358115051988285458967698823);
[ <634466267339108669, 1>, <3865430919824322067, 1> ]
> primes := [ f[1] : f in Factorization(p-1) ];
> p1, p2, p3, p4, p5, p6 := Explode(primes);
> exponents := [ (p-1) div r : r in primes ];
> n1, n2, n3, n4, n5, n6 := Explode(exponents);
> FF := FiniteField(p);
> a := FF!109;
> c := FF!1014452131230551128319928312434869768346;
\end{verbatim}
